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Abstract 


The bimaximal (BM) neutrino mixing matrix was formulated in order to accommodate the data of 
the experimental results which indicate that both solar and atmospheric neutrino oscillation in vacuum 
are near maximal. But, after the T 2 K and Daya Bay Collaborations reported that the mixing angle 613 is 
nonzero and relatively large, many authors have modified the neutrino mixing matrix in order to accom¬ 
modate experimental data. We modified the BM mixing matrix by introducing a simple perturbation 
matrix into BM mixing matrix. The modified BM mixing matrix can proceed the mixing angles which are 
compatible with the globat fit analysis data and by imposing the A* ~symmetry into mass matrix from 
modified BM, we have the neutrino mass in normal hierarchy: mi < m2 < m3. Using the neutrino masses 
that obtained from neutrino mass matrix in the scheme of modihed BM and imposing the constraint 
exact — r symmetry into neutrino mass matrix, we cannot have compatible squared-mass differences 
for both Am|i and Am32 as dictated by experimental results. We break softly the M — t symmetry by 
introducing a small parameter A into neutrino mass matrix which then can proceed neutrino masses are 
in agreement with the squared mass difference as dictated by experimental results.The predicted neutrino 
effective mass: |meel = 0.0155 eV in this paper can be tested in the future neutrinoless double beta decay. 
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1 Introduction 

Neutrino oscillation phenomena can be explained if neutrino has a nonzero mass and some mixing angles 
exist in neutrino sector. The mixing angles for three neutrinos can be formulated in mixing matrix. The 
mixing matrix relate the neutrino eigenstates in flavor basis {u„, into neutrino eigenstate in mass 

basis {mi,m2,m3) as follow: 



( 1 ) 


where V is the mixing matrix. The standard parameterization of the mixing matrix read: 



where dj is the cos 6ij, Sij is the sinSij, 6ij are the mixing angles, and 5 is the Dirac CP-violating phase. 

There are three well-known of neutrino mixing matrix, i.e. bimaximal mixing (BM), tribimaximal 
mixing (TBM), and democratic mixing (DM). All of the mixing matrices predict the mixing angle 613 = 0 . 
Recently, the experimental results showed that the mixing angle 613 yf 0 and relatively large Bill El mg. 
Several attempt have been done theoretically to accommodate the nonzero and relatively large mixing 
angle O13 by modifying the neutrino mixing matrix including the Dirac phase S in relation to the CP- 
violation in neutrino sector. Another unsolved problem in neutrino physics till today is the hierarchy of 
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neutrino mass. Experimental results showed that we have two possibilities for neutrino mass hierarchies: 
normal and inverted hierarchies. We have no clue in order to decide theoretically the neutrino mass 
hierarchy whether it normal or inverted. 

Theoretically, the neutrino masses can be determined from neutrino mass matrix (Mi,). The neutrino 
mass matrix is related to the neutrino mixing matrix via the following equation: 

M„ = EME^, (3) 

where M is the neutrino mass matrix in mass basis and V is the mixing matrix. From Eq. , we can see 
that the neutrino mass pattern in flavor basis depend on the pattern of the neutrino mixing matrix. We 
also belief that the pattern of the neutrino mass matrix should be related to one of the unique underlying 
symmetry. One of the interesting candidate symmetry is the /r — r symmetry because it can reduce the 
number of parameters in neutrino mass matrix and can predicts qualitatively the neutrino mass hierarchy 
that compatible with the experimental results. But, if we use the advantage of the experimental results as 
input to fix the values of some parameters that we build from theoretical side with exact /r — r symmetry, 
many theoretical predictions are incompatible with the experimental results especially mixing angle ^13 
which lead to be zero in the frame of ^ — r symmetry. Thus, if we still want to use the fJ. — r symmetry, 
the we should invoke a perturbation into neutrino mass matrix or into three well-known mixing matrices 
BM, TBM, and DC. The idea of introducing a perturbation into — t symmetry mass matrix have 
been introduced in la [71 n |9], where authors analyzed the effect of perturbation and its correlation 
corresponding to the mixing angles: O 12 and 623 . In m the broken fj, — t symmetry is used to obtain 
nonzero 613 and Jarlskog rephasing invariant. 

In this paper, we determine the neutrino mass hierarchy and neutrino masses from the neutrino mass 
matrix that obtained from the modified neutrino mixing BM with additional constraint softly broken 
/i — r symmetry as the underlying symmetry of the resulted neutrino mass matrix. In section II, we 
modified BM by introducing a simple perturbation matrix and calculate mixing angles 6^12 and 623 by 
taking the advantage of Daya Bay Collaboration result on mixing angle 613 H]. In section III, we evaluate 
the neutrino mass hierarchy and neutrino masses from the neutrino mass matrix obtained from modified 
BM with assumption that the underlying symmetry of the neutrino mass matrix is the fJ. — r symmetry. 
We also discuss predictions of the obtained neutrino masses especially on the squared mass difference for 
both solar and atmospheric neutrinos. Finally, section IV is devoted for conclusions. 


2 Modified BM Mixing Matrix 

The bimaximal (BM) mixing matrix was formulated in order to accommodate the facts that both solar 
and atmospheric data can be described by maximal mixing vacuum oscillation with the relevant mass 
scale and it imply that there is a unique mixing matrix which is then called bimaximal mixing matrix. 
The BM mixing matrix {Vbm ) read |11 |: 


/ \/7 0 \ 

VsM = I 5 5 -\/T 1 ■ (^) 

Vi i ViJ 

As one can see from Eq. @, by comparing it to neutrino mixing matrix in Eq. m, the BM mixing 
matrix give: sin^ise”*"^ = 0 which imply that the BM mixing matrix leads to mixing angle 833 = 0 which 
is incompatible with the recent experimental results of T2K Collaboration pQ: 

5° < 6 »i 3 < 16°, (5) 

for neutrino mass in normal hierarchy (NH), and 

5.8° < 6 »i 3 < 17.8°, ( 6 ) 

for inverted hierarchy (IH) with Dirac phase: <5 = 0. The nonzero value of mixing angle ^13 was also 
conhrmed by Daya Bay Collaboration as follow [4]: 

sin^ 26 »i 3 = 0.092 ± 0.016 (stat.) ± 0.005 (syst). (7) 

From Eq. we also see that the mixing angle: 0 i 2 = 823 = 7 r /4 (bimaximal mixing). 

In order to accommodate the relatively large and nonzero mixing angle 813 in accordance with BM 
mixing matrix, many authors have modihed the BM mixing matrix |1211131114 | . In this paper we modify 
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the neutrino matrix mixing BM by introducing a simple perturbation matrix to perturb BM which 
different with the Refs. |12II131[T^ . The perturbation matrix {Vp) is given by (with Dirac phase: <5 = 0): 


Vp = 



0 

1 

0 



( 8 ) 


where Cp is cosp, and Sp is sinp. The modified BM mixing matrix {Ubm) is obtained via the following 
relation: 


which then gives: 


Ubm — VpVsM, 


Ubm = 



(9) 


( 10 ) 


From Eq. GQl) one can see that the modified BM can predicts nonzero mixing angle 613 . By comparing 
Eq. (nni) with Eq. ©, we have: 


tan 012 = 


Sp ^^^Cp 
Sp + \picp 


tan 023 = 



sin 013 = 


^\f^Sp 

2 


( 11 ) 

( 12 ) 

(13) 


By inspecting Eqs. (HU, (HSJ, and m, we can see that the mixing angles: 0i3 yf 0, 0i2 < 7r/4, and 
023 >- 71 / 4 . On the other hand, the modified BM in this scenario proceed no-maximal mixing matrix. 

If we insert the central value of mixing angle 0i3 as reported by Daya Bay Collaboration as written 
in Eq. ©, then we have: p = 12.54° and it implies mixing angles: 

012 = 36.07° and 023 = 45.69°, (14) 

which is compatible with the global fit analysis data |15 |: 

012 = 34.5 ±1.0(ll j)°, (15) 

023 = 42.8l®;®(t^°3")°, (16) 


for Icr {3a) level. 


3 Neutrino Mass Hierarchy from Ubm 

As stated in section I, one of the the unsolved problem in neutrino physics is the neutrino mass hierarchy 
whether it normal or inverted. By using the Eq. with neutrino mixing is the modified BM {Ubm) 
with mixing angles as shown in Eq. (HI and mass matrix (M) in mass basis is diagonal: 


/ mi 0 0 \ 

M = j 0 m 2 0 j, 
Vo 0 msj 


then we have the neutrino mass matrix as follow: 


3 follow: 




Q 

R 

= Q 

s 

T 

\R 

T 

W 


(17) 


.3 


M. 


(18) 














where: 


P = 0.638mi + 0.338m2 + 0.024m3, (19) 

Q = —0.904mi + 0.084m2 + 0.108m3, (20) 

R = -0.115ml - 0.5764m2 + 0.106m3, (21) 

S = 0.979mi + 0 . 021 m 2 + 0 . 499 m 3 , (22) 

T = —0.143mi + 0.143m2 + 0.488m3, (23) 

W = 0.021mi + 0 . 979 m 2 + 0 . 477 m 3 . (24) 


If we impose the exact — t symmetry into resulted neutrino mass matrix (M^) of Eq. (18), then 
we must put Q = R or 


and S = IE or 


Q- R = 0, 

S' - IE = 0. 


Solving simultaneously Eqs. (I25l) and (1261) . we have neutrino masses relations as follow: 

1712 ~ 1.1727mi and m 3 = 7.5190mi, 


(25) 

(26) 

(27) 


which imply that: 

mi < m 2 < m 3 . (28) 

The resulted hierarchy of neutrino mass in Eq. (1281) is normal hierarchy (NH). 

By taking the advantage of global fit analysis of the experimental results of the squared-mass difference 

Amh = 7.59 ± 0.20(to:69) x 10"® eV^ (29) 

Am 32 = 2.46 ± 0.12(±0.37) x 10”® eV^ (for NH), (30) 

then we have: 

mi = 0.014224 eV, 

m 2 = 0.016680 eV, (31) 

m 3 = 0.106949 eV, 

when we use the squared-mass difference of Eq. (1291) to determine the neutrino mass in Eq. dSTJ. 
Conversely, if we use the squared-mass difference of Eq. (1301) to determine the neutrino mass in Eq. (1271) . 
then we have: 


mi = 0.006678 eV, 

m 2 = 0.007831 eV, (32) 

m 3 = 0.050212 eV. 

From Eq. ®I), one can see that the neutrino masses predict the squared-mass difference Am ^2 ~ 
11.16 X 10“® eV^ which is incompatible with the experimental result. It is also apparent from Eq. ([32]) 
that neutrino masses predict Am|i = 1.673 x 10“® eV^ which is incompatible with the experimental 
results. Thus, in the context of exact — t symmetry as an underlying symmetry of neutrino mass 
matrix, neutrino masses cannot proceed compatible predictions for both squared-mass differences:Am 2 i 
and Am §2 as dictated by experimental results. In order to get the compatible predictions of the obtained 
neutrino masses with the experimental results, we break softly the fJ. — t symmetry with the following 
scenario: Q = R or 


Q- R = 0, 


(33) 


and S - IE = A or 


S' - IE - A = 0, (34) 

where A is a small parameter that perturb softly the exact — t symmetry. Solving Eqs. ([331) and m 
simultaneously, we have: 


mi = 0.1329967m3 - 5.3405836A, (35) 

m 2 = 0.1559612m3 - 6.3844249A. (36) 


A 


Using the neutrino mass relations in Eqs. (|35ll and (I36II and we equate it with the squared-mass 
difference of Eq. (I29II to determine the neutrino mass m 3 , then we have: 

m 3 = 43.0157073A -f 0.964468774 x 10"^®\/0.63942498A2 -f 0.122962768 x lO^s. (37) 

Inserting the global fit analysis result in Eq. (30) into squared-mass difference Am §2 which is calculated 
from Eqs. IMl) and ([37ll. we have the value of parameter A as follow: 

A = -0.00127381 or A = -0.00370340. (38) 

By inserting the value of A = —0.00127361 into Eqs. (I35II . (I36II . and (I37II . we finally have neutrino 
masses: 


mi = 0.0137452 eV, 

m 2 = 0.0162737 eV, (39) 

m 3 = 0.0521999 eV. 

that can proceed the squared-mass differences Am|i and Am 32 are compatible with the experimental 
results. By using the relation of neutrino effective mass of neurinoless double beta decay: 

\mee\ = \ {UBM)eimi + (?7flM)e2”1-2 + (UflM)e3m3| , (40) 

then we have: 

|meel = 0.0155 eV, (41) 

that can be tested in the future neutrinoless double beta decay experiment. 


4 Conclusions 

We have modified BM mixing matrix by introducing a simple perturbation matrix into BM that can 
proceed nonzero mixing angle: 613 and relatively large. From the modified BM mixing matrix and using 
the central value of mixing angle: ^13 from Daya Bay Collaboration data, we determine the mixing angles: 
012 and 023 which are compatible with the global fit analysis data. The neutrino mass matrix obtained 
from the modified BM mixing matrix with is constrained by exact symmetry proceed the neutrino mass 
in normal hierarchy: mi < m 2 < m 3 , but neutrino masses cannot predict the squared-mass differences: 
Am^i and Am 32 as dictated bay the experimental results. To get the neutrino masses that can proceed 
the compatible squared-mass differences: Am|i and Am 32 with the experimental results, we break softly 
the — T symmetry by introducing a small parameter A into neutrino mass matrix. The obtained 
neutrino masses from softly broken fJ, — T symmetry as the underlying symmetry of neutrino mass matrix 
can proceed the squared-mass differences: Am|i and Am §2 which are in agreement with the experimental 
results with neutrino masses: mi = 0.0137452 eV, m 2 = 0.0162737 eV, and m 3 = 0.0521999 eV. The 
predicted neutrino effective mass: |meel = 0.0155 eV in this paper can be tested in the future neutrinoless 
double beta decay 
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